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We study the simplest Two-Higgs-Doublet Model that allows for CP nonconser- 
vation, where it can be parametrized by only one parameter in the Higgs poten- 
tial. Different concepts of maximal CP-nonconservation in the gauge-Higgs and the 
quark-Higgs (Yukawa) sectors are compared. Maximal CP nonconservation in the 
gauge-Higgs sector does normally not lead to maximal CP nonconservation in the 
Yukawa sector, and vice versa. 

1 Introduction 

Mendez and Pomarol introduced the concept of maximal CP nonconservation [1] in the 
context of the gauge-Higgs sector of the Two-Higgs-Doublet Model (2HDM) [2]. In the 
absence of CP nonconservation, only two of the three neutral Higgs bosons couple to the 
electroweak gauge bosons (the two CP even ones, often denoted h and H). When CP is not 
conserved, all three do. In fact, Mendez and Pomarol realized that the product of all three 
gauge-Higgs couplings, which is bounded by unitarity, is a useful concept to parametrize 
the amount of CP nonconservation, and defined the quantity 
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£y — 27 [gvv h x 9vvh 2 9vvh 3 ] 
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as a measure of CP nonconservation in the gauge-Higgs sector. If the couplings gvvHi are 
normalized with respect to those of the Standard Model (SM), then £y, as defined above, 
satisfies 

< & < 1. (1-2) 

However, this measure of CP nonconservation is not applicable to the fermion-Higgs sector. 

In the fermion-Higgs sector of a given version of the 2HDM, one should consider quanti- 
ties other than £y as measures of CP nonconservation. As we will see from our investigation, 
the parameters of the 2HDM that maximize £y are different from those that maximize CP 
nonconservation in the Yukawa sector. They are in general also different for the up- and 
down-quark sectors. 

The paper is organized as follows. In sect. 2 we review the 2HDM and in sect. 3 we 
study the conditions for maximum CP nonconservation in the gauge-Higgs sector. Sections 
4 and 5 are devoted to the Yukawa sector, at the parton and proton level, respectively, and 
sect. 6 contains some concluding remarks. 



2 The Two-Higgs-Doublet Model 

We shall here introduce some notation for the Two-Higgs-Doublet Model [3]. Let the Higgs 
potential be parametrized as [4] 



V = y(0 t 1 0i) 2 + y(40 2 ) 2 + A 3 (0 t 1 0i)(0 2 2 ) + A 4 (0 t 1 2 )(40 1 ) 

1 r 



(2.1) 



+ 2 



A 5 (4 



+ h.c. 



+ 



m i 2 (0i02) +h.c. 



+ m 2 22 



(402)}. 



The parameters A 5 and m\ 2 are allowed to be complex, subject to the constraint 

Im TO12 = ImA 5 i>ii>2, (2.2) 

with v\ and v 2 the vacuum expectation values (vf + v\ = v 2 , with v = 246 GeV). 
The corresponding neutral-Higgs mass matrix squared is then given by 

/ Aic^ + z/s^, (A 345 - v)cpsp -\lm.\ b sp\ 

M = v 2 (A345 -v)cpsp \ 2 s 2 p + vc 2 p -ilmAsC^ (2.3) 

\ — |lmA 5 s^ — |lmA 5 c^ — Re A 5 + v) 
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with the abbreviations cp = cos/5, sp = sin/3, tan/5 = v^/vi, A345 = A3 + A4 + ReA5, 
v = Rem^ 2 /(2v 2 sin/5cos/9) and /i 2 = v 2 u. 

The (1,3) and (2,3) elements of this mass-squared matrix (2.3), which are responsible 
for CP nonconservation, are related via the angle (3. In this sense, CP nonconservation is 
described by one parameter, namely ImA 5 . 

In order to diagonalize this matrix (2.3), we introduce the rotation matrix 





R — R c R b R & = I cos a c 




( cos a b sina fe N 

1 
\— sin a b cos a b/ 



( cos a sin a s 
- sin a cos a 
1, 



-(c & s b s c + S 5 C C ) CaC c -S & S b S c c b s c \ (2.4) 
Sb c c + s b c c ) C b C C/ 



with Cj = cosaj, Sj = sinctj, and satisfying 

RM.R T = diag(M 1 2 , M 2 , M 2 ). (2.5) 

Here, M 1 < M 2 < M 3 . The angular ranges are taken as — n/2 < a < ir/2, —ir < a b < 7r, 
and — 7r/2 < a c < 7r/2. As discussed in [5], only some regions of the parameter space are 
physically allowed. 

This limitation of the parameter space is due to various constraints, including (i) Mi < 
M 2 < M 3 , and (ii) the constraints of perturbativity and unitarity. We shall represent the 
latter as 

I Aj| < 47r£pcrt, With £pcrt = 0(1). (2.6) 

We show in Fig. 1 typical allowed regions in the a b -a c plane, for a few values of tan f3 
and a. In this figure, we only show regions of |q;j,| < 7r /2 and only positive a c . Regions 
of larger |cKbj and negative a c are given by the symmetries discussed in [5]. Furthermore, 
for given values of tan/5 and a (and given sign of a c > 0), only one sign of a b is realized, 
requiring M 2 < M 3 . The dashed lines at a b = ±7r/4 indicate where CP nonconservation 
is maximal in the Higgs-top-quark sector, in the limit of one light Higgs boson and two 
heavier ones, see Eq. (5.4). 

Different choices for the 'soft parameters' (in particular, different values of fi 2 ) lead 
to somewhat different allowed regions. Also, a larger value of £ per t extends the region. 
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Figure 1: Dark (blue): physical regions [see Eq. (2.6)] in the ab-a c plane for various values 
of tan (3 and a. Soft parameters: M x = 100 GeV, M 2 = 500 GeV, M H ± = 600 GeV, 
ji = 300 GeV, £ pcrt = 1. Light (yellow): Same with £ pert = 5. Solid contour: absolute 
boundary. 



However, there are absolute bounds, indicated by the solid contours outside the shaded 
regions in Fig. 1, that can not be crossed for any choice of the 'soft parameters' [5]. In 
order to cover a range of different choices for /i 2 , one may take a rather large value of £ pcrt 
(in sect. 5 we shall consider £ pert = 5). For further discussion of these issues, see [5,6]. 

In this notation, Eqs. (2.3)-(2.5), the gauge-Higgs couplings are, relative to the corre- 
sponding SM coupling, given by 

KiZZ : 9vVHi = cos(3Rn + sin/3i? i2 , (2.7) 
whereas for the Yukawa couplings we consider the so-called Model II [3] where they are 
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given by 

Hjti : -— [R j2 - h 5 cos PR j3 \ = [af + i-y b af\, (2.8) 
sin 13 

H^bb : — *— [% - 2 75 sin /^ 3 ] = [af + i 75 af ] , (2.9) 

COS /./ 

with i?^ an element of the rotation matrix (2.4). 

3 CP nonconservation in the gauge-Higgs sector 

In the gauge-Higgs sector, the amount of CP nonconservation [cf. Eq. (1.1)] is in the above 
notation given by 

3 

£v = 27 Y[[cosPRn + sm(3R t2 } 2 . (3.1) 
i=i 

This £y depends on tan (3 as well as on the three angles a, aj, and a c that determine Rij. 
However, it only depends on (5 and a through their difference. In fact, using (2.4) and 
some trigonometric identities, we find 

£y = 27c^cos 2 (/9 — a)[sbS c cos((3 — a)— c c sin(/9 — a)] 2 [sbC c cos(j3 — a) + s c sm({3 — a)] 2 . (3.2) 

It is also seen that £y is unchanged under 

a b fixed, (a c <-> 7r/2 + a c ) : £y <-> £y, 

(a 6 <-> (a c <-> 7r/2 - a c ) : £y <-> £y. (3.3) 

In order to provide some intuition for how the CP nonconservation depends on the 
parameters of the 2HDM, we show in Fig. 2 contours of constant £y in the ai,-a c plane, for 
various values of tan/3 and a. We note that there is little CP nonconservation for 'large' 
values of because of the factor (? h in (3.2). Also, there is CP nonconservation even for 
«6 = and for a c = (but not when both vanish). 

3.1 Simple limits 

It is instructive to consider the simple limits of a b = or a c = 0. 
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tan/3 = 0.5 tan/?=1.0 tan/? = 2.0 




Figure 2: Contours of constant £ v [see Eq. (3.1)] in the ab-a c plane for various values 
of tan p and a. Soft parameters: M l = 100 GeV, M 2 = 500 GeV, M H ± = 600 GeV, 
H = 300 GeV. Dark (blue): £ pert = 1, light (yellow): £ pcrt = 5. 



Oi b = 

For otb = 0, the rotation matrix simplifies: 



R = I -Sa c c c & c c s c I , (3.4) 
and one finds 

27 

£v{(*b = 0) = — sin 2 (2a c ) sin 4 (/5 — a) cos 2 (/5 — a). (3.5) 
The maximum is given by 

£y — 1 for a = /5 ± arctan v^2, = 0, a c = ±7r/4. (3.6) 
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For a c = 0, one finds 

27 

£y(ai c = 0) = — sin 2 (2a;,) cos 4 (/3 — a) sin 2 (/3 — a). (3.7) 
This relation holds also for a c = ir/2. The maximum is given by 

£y = 1 for a = /3 ± arctan(l/v / 2), = ±7r/4, a c = or a c = tt/2. (3.8) 
3.2 Maxima of £ v 

Since maximizing over angles allows us to keep two Higgs masses fixed [5] and since by 
Eq. (3.2), the dependence of t;v on P an d ol shows up in the form [j3 — a), £y can be 
maximized for fixed {(3 — a) by meeting the two conditions: 

|K = and fK = . (3.9) 

By substituting from Eq. (3.2), and solving (3.9) for ay, and a c , we obtain a continuum of 
maxima: 



l+tan 2 (/j-q) 
4y = 1 tor a;, = ± arccos y , 

1 + tan 2 (/3 - a) - a/3[2 - tan 2 (/5 - a)l tan(/3 - a) . 

a c = ± arctan 0//n — . 3.10 

2tan 2 (/3- a) - 1 v ; 

which impose the constraint 

\ttm((3-a)\ < V2 (3.11) 

on {j3 — a). We note that (3.6) and (3.8) are both special cases of this (3.10). 

We show in Fig. 3 how these angles ab and a c vary with tan (5 (for fixed a) when we 
maximize £y. For a given value of <S, these curves only cover a finite range in tan/5. They 
are cut off by (3.11), which says that, in order to have £y = 1, f3 and a should not differ 
by more than arctan \[2 ~ 54.7°. In addition, they are cut off by the condition of having 
a physical solution as discussed in sect. 2, and delineated by the solid contours in Fig. 1. 
Note that there are also solutions having other signs for ct b and a c , but that the model is 
only physically consistent for certain sign combinations. 
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Figure 3: Angles ab and a c [cf. Eq. (3.10)] for which the CP nonconservation £y in the 
gauge-Higgs sector is maximal, for a range of tan/? values, and for a = 0, tt/6, 7r/4, and 
?r/3. 



4 CP nonconservation in the Yukawa sector 

In the Yukawa sector, one can define measures of CP nonconservation analogous to the 
one for the gauge-Higgs sector [cf. £y of Eq. (1.1)]. Requiring thus that all three Higgs 
bosons should have CP-nonconserving couplings to up- and down-type quarks, it is natural 

to consider the quantities [see Eqs. (2.8), (2.9) and (5.3)]: 

6 3 / ___ a x 6 

7t, 



\sm p J rr \ sin p J 



Both of these differ from the defined above in two respects. First of all, the dependence 
on (3 factorizes. Secondly, they individually diverge as sin (3 — > (for up-type quarks) or 
cos 13 — > (for down- type quarks) . 

One could also consider the quantities 

sin p J r-f 

' ' i=i 

and similarly as measures of CP nonconservation in the Yukawa sector. These measures — 
unlike those in (4.1) — are consistent with the fact that if H\ conserves CP in its couplings 

8 



to the up- and down-type quarks, i.e. a b = 0, then the Yukawa sector may still be CP 
nonconserving, since the other two Higgs states, H 2 and H 3 , may have CP nonconserving 
couplings to the quarks. Accordingly, ( t 7^ and C& 7^ f° r a b = which is not the case 
for £ t and £&. This ( t will be discussed in sect. 4.3. 

Substituting now from (2.4) into (4.1), we obtain for this case of Model II Yukawa 
couplings: 

Tt ^b i^a $b Cc S c ) [s<j C c + Cq S c ] [s<j S5 S c C5, C c ] , 

7b = Cfe(c 5 s fe c c s c ) 2 [c fi s b s c + s & c c ] 2 [c & s b c c - s & s c } 2 . (4.3) 

We note that both these quantities possess the same symmetries (3.3) as £y. Also, % is 
obtained from j t by the substitutions 

(sa <-»■ Ca), (s c <-> Cc) : 7t <-> %■ (4.4) 



4.1 Maxima of j t 

Let us now consider the maxima of % in (4-3). We find the maximum value ^™ ax = 1/1024 
for 

Case I : a — ^n, a b = ±\n, a c = ±\n, (4.5) 
where the two signs are independent, and at 

Case II : a — ±arctan^= (at = ±0.196 7r), a b = ±|vr, with 

V2 

a c = iarctan^p (a c = ±0.196 ir) or a c = ^arctanv^ (a c = ±0.304 ir). (4.6) 
v2 

For Case II, the signs are subject to the constraint aa b a c > for the first a c solution, and 
aa b a c < for the second a c solution. The maxima of % are obtained by the substitutions 
(4.4). 

Thus, it is natural to define normalized quantities 

3 3 

lt = 1024 ]J[R i2 R i3 ] 2 , lb = 1024jJ[i2 il R l3 ] 2 , (4.7) 
i=\ i=i 
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Figure 4: Contours of constant 7 t [see Eq. (4.7)] in the a&-a; c plane for various values 
of tan p and a. Soft parameters: M 1 = 100 GeV, M 2 = 500 GeV, M H ± = 600 GeV, 
fj, = 300 GeV. Dark (blue): £ pert = 1, light (yellow): £ pcrt = 5. 



satisfying 



< 7t < 1, < 7 6 < 1, 



as measures of CP nonconservation in the up- and down-quark sectors, respectively. Con- 
tours of constant j t are shown in the afc-a c -plane in Fig. 4. 
Let us now keep a fixed. Then, the maxima of jt are at 



Case I : a& = e b |7r, a c = e c arctan v^-^tan 2 a + \ + e&e c tan 1 
Case II : cn&= e;, |7r, a c = e c arctan 



a 



| ( a/ tan 2 a + A — e b e c tan a 



(4.9) 



where e& and e c are independent sign factors: = ±1, e c = ±1. For Case I, the corre- 
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sponding maximum is (same for all sign choices) 

7t = sin 6 a, (4.10) 

in agreement with Eq. (4.5), whereas for Case II, the corresponding maximum is (same for 
all sign choices) 

27 tan 2 a , . 

7t = T (l + tan 2 «)3' (4 - H) 

which becomes 1 for tana = ±l/\/2, in agreement with Eq. (4.6). 
4.2 Maxima of £ Y 

While £ t and individually diverge as (3 — > and (3 — > 7r/2, respectively, the product over 
couplings to up-type and down-type quarks is less divergent. We define, analogous to (1.1) 
and (4.1) 

(412) 

with 

3 

1y = 7o It lb = 7o J] [ R * R ^ R U 2 ■ ( 4 13 ) 

1=1 

satisfying 

< 1y < 1. (4.14) 

Substituting from (4.3), we obtain 

7y = 7o c^ 2 (c c s c Sfe) 4 (c 5 s<i) 2 [s<5 c c s fe + c 5 s c ] 2 [os c c s 6 - s fi s c ] 2 

x [c 5 s b s c + Sa c c ] 2 [s fi s 6 s c - c 5 c c ] 2 . (4.15) 

This has a maximum for (see Appendix A) 

S = ±±7T, a 6 = ±arcsiny // | = ±0.13386tt (24.1°), a c = ±\n, (4.16) 

with 



^ = ^Ttr-=( ^ J = 3.652 xlO 6 . (4.17) 



2 26 3 i2 _ /8 x 1024 x 27 2 \ 2 
5 10 ~ V 3125 

Fig. 5 exhibits contours of constant 7y for some values of a other than that of the 
maximum, a = ^n, in relation to the physically allowed (dark, shaded) regions in the 

11 



tan/3 = 0.5 tanjff=1 .0 tan/? = 2.0 




Figure 5: Contours of constant 7y in the a(,-a c plane for various values of tan/5 and a. 
Soft parameters: M x = 100 GeV, M 2 = 500 GeV, M H ± = 600 GeV, fi = 300 GeV. Dark 
(blue): ^pe r t = 1, light (yellow): £ pcTt = 5. 



atb-a c plane. Note that 7y vanishes when a = or a = ±7r/2, as well as on the edges of 
the quadrants: = or ±7r/2, a c = or ±7r/2. Also, we note that there are secondary, 
local, maxima. 

Although 7y is, by definition, independent of (3, Fig. 5 shows contours of constant 
7y superimposed on allowed regions for different values of tan/?, since the 'shapes' and 
locations of the physically allowed regions in the ol\>-oi c plane depend on tan f3. Accordingly, 
the positions of the maxima 2 of jy, w.r.t. the physically allowed regions in the af,-a c plane 
are different for different values of tan/3. For example, consider a = tt/Q. We see from 
2 This is not a 'maximum' in the same sense as above, since a is held fixed. 
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Fig. 5 that for tan/3 = 0.5, 7™ ax is located outside the physically allowed region while for 
tan (3 = 1.0, this is not the case. Moreover, for tan/5 = 2.0, the physically allowed region 
shifts the location to the 'other' quadrant. To sum up, for a — 7r/6, the location of 7™ ax 
occurs at 

(a b , « c )| tan /3=0.5 = («&, «c) | tan 0=1.0 = (-a&,7r/2 - OL c ) | tan 0=2.0- 

4.3 Maximizing Ct 

We now return to the quantity Ct of Eq. (4.2), which we rewrite as 

C t = (cos [3/ sin 2 [3) 2 ( t (4.18) 

with 

3 

Q = 2^[ J R j2 i? l3 ] 2 , < Ct < 1- (4-19) 

i=i 

Substituting from (2.4) and utilising trigonometric identities, we find 

Ct = \c 2 h [{l - C 2& )(7 + C Ac )s 2 h + 2s 2& S 4c S b + (1 + c 2& ){\ - c 4c )]. (4.20) 

To maximize Ct, we differentiate w.r.t. a, a b and a c and get: 

S2a(7 + C Ac )s 2 b + 2c 2(5 S4 C Sb - «2a(l ~ C 4c ) = 0, 

2(1 - c 25! )(7 + c 4c )s^ + 3s 2& s 4c s 2 b - 2(3 - Ac 2& + c 4c )s b - s 2& s 4c = 0, 

(1 - c 2& )s 4c s 2 b - 2s 2& c 4c s b - (1 + c 2fi )s 4c = 0. (4.21) 

Solving the three equations, one finds: Ct — 1 f° r 

Case I : c 2 a = 1, s b — 0, c 4c = —1 

Case II: c 25l = -1, s 6 = il/v 7 ^, c 4c = 1 (4.22) 

with the corresponding angles 

Case I : a — 0, = or oj = ±7r, a c = ±|7r, 
Case II : a — ±§7r, a: 6 = ±|7r, a c = or ± |/T. (4.23) 
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tan/3 = 0.5 tanjff=1 .0 tan/3=2.0 




Figure 6: Contours of constant Q [see Eq. (4.20)] in the af,-a c plane for various values 
of tan p and a. Soft parameters: M 1 = 100 GeV, M 2 = 500 GeV, M H ± = 600 GeV, 
fj, = 300 GeV. Dark (blue): £ pert = 1, light (yellow): £ pcrt = 5. 



Considering now a fixed, we find the maxima: 

Case I : a.\, = (or ± it), a c = ±t7T, 

for which 

Ct = |(l + C2a) 

coincides with Case I in (4.22) for c 2 a = 1, and 

V3 - 5c 2 a 



Case II : = ± arcsin 



a, = ±- 



7r — arccos ■ 



5 + 13c 2 a 
11 + 3c 2 5 



(4.24) 



(4.25) 



(4.26) 
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provided C25 < 3/5. In this case 

6 = 32(1 -l^^ K" + 3Cm)(43 - 10C - + U 4> 
+ (44 + 12c M ) V3 - 5c 2 «VTTS\/3 - 2c M - 54J , (4.27) 

agrees with Case II in (4.22) for C25 = —1. 

Fig. 6 exhibits contours of constant (t in the a b -a c plane for selected values of tan/3 
and a. We read off from Fig. 6 that for a — 0, the quantity ( t takes its maximum value 
at (a b ,a c ) = (0, \n) which again is consistent with Case I in (4.22). For particular values 
of a and a c , there are also saddle points, for example at (a, a b , a c ) = (|-7r, — j7c, jn). For 
a given value of a, these saddle points are located at 

a b = ±±7r, a c = ±\ arccos ( + C2a ) (4.28) 

V 3 + c 2 & J 

in the a b -oi c plane. 

In the top-Higgs Yukawa sector, 7 t [see Eq. (4.7)] and ( t are both sizable for large a and 
\a b \ ~ vr/4, as we see in Figs. 4 and 6. However, the two measures have different features. 
For example, for the same value of a — 0, where ( t has a maximum, 74 vanishes. Moreover, 
for a b = 0, 7 t vanishes (regardless the values of <5 and a c ) while ^ takes its maximum value 
(for o: = and a c = This again shows that these quantities 74 and behave rather 
differently for a given set of the angles (a, a b , a c ). 

5 CP nonconservation in pp — > ti 

The above studies refer to the tree-level couplings of Higgs particles to vector particles 
and fermions. These are difficult to study directly, since the Higgs particles as well as 
the vector particles and the relevant fermions are unstable. The implication is that it is 
easier to access these couplings via various loop effects. We shall here consider one such 
example, namely the production amplitudes for the tt through gluon fusion, where CP 
nonconservation is induced by non-standard neutral Higgs exchange. 

CP nonconservation in the production of tt pairs at future hadronic colliders has been 
studied in considerable detail [7]. For a detailed application to the 2HDM, see also [5]. 
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One process of particular interest is 



pp^ttX, (5.1) 

where the t and t decay semileptonically, and the lepton energy difference is measured [5,7]: 

Ai — E + — E_. (5.2) 

(For a discussion of other observables, see [7,8].) The expectation value of this observable 
will in general be non-zero if there between the quarks in the final state are exchanges of 
Higgs bosons that are not eigenstates under CP. The quantity [see Eq. (2.8)] 

lcp,j - -a, a j - — ^ K j2 K j3 (5.3) 

then plays a crucial role, together with non-trivial functions of the kinematics (given by 
the loop integrals). 

If the neutral-Higgs spectrum has a large gap between the lightest Higgs boson and the 
next one, then the lightest one will give the dominant contribution to Ai, and the amount 
of CP nonconservation is roughly proportional to 



1 sin a sin(2afe) 
tan (3 sin (5 



7CP.1 = 2 — a > ( 5 - 4 ) 



which is maximized for small tan/5 and for (<5, a&) = (±7r/2, ±7r/4), corresponding to 
the dashed lines at = ±7r/4 in Figs. 1, 2, 4, 5 and 6. These values as well, (a, c*&) = 
(±7r/2, ±7r/4), together with a c = or ±§7r, coincide with those of Case II [see Eq. (4.23)] 
that maximize ( t . Furthermore, (a, a^) = (tt/2, ±7r/4), together with a c = ±|7r, coincide 
with those of Case I [see Eq. (4.5)] that maximize 7 t . These results indicate that large 
o: together with ~ 7r/4 favour large CP nonconservation in the Yukawa sector. It 
is immediately obvious that this is not compatible with the condition of maximal CP 
nonconservation in the gauge-Higgs sector [1], £y = 1 [see Eqs. (3.6) and (3.8)]. 

In addition to the contribution from the lightest Higgs boson, there will in general 
also be non-negligible contributions from the others. Because of the orthogonality of the 
rotation matrix R, not all Jcpj can have the same sign, so there will be cancellations. 
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Let us define the 'signal-to- noise ratio', or sensitivity [7] 

S = (A,) 

N 

which provides a measure of how much data would be required to see an effect. 



(5.5) 
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Figure 7: Left panel: Maximal sensitivity [see (5.5)] for the observable (5.2), for fixed M±, 
M 2 and two values of tan (3. Right panel: Corresponding values of the angles a, and a c . 
Soft parameters: M 2 = 500 GeV, M H ± = 600 GeV, £p Cr t = 5. 



It is interesting to maximize the amount of CP nonconservation that results for the 
observable Ax, over the relevant parameters of the model. In Fig. 7 we show the result 
of such a maximization of the sensitivity (5.5). The quantity A± and its spread A\ are 
computed as given in [5,7], using the 'LoopTools' package [9, 10], and convoluted with the 
CTEQ6 parton distribution functions [11] for the LHC energy of 14 TeV. The resulting 
quantity is then maximized using the 'MINUIT' package [12]. 

The actual maximization is rather CPU-intensive: In order to evaluate A\ and S/N, 
three-dimensional integrals (a convolution integral over the parton distribution functions, 
an integral over the polar angle of the top quark with respect to the beam, an integral 
over s, the invariant mass squared of the ti pair) involving non-trivial loop functions are 
required. These are then maximized in the three angles parameterizing the 2HDM mass 
matrix: a, a 6 and a c (keeping the two lowest Higgs masses fixed). 
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In this maximization, we have kept M 2 = 500 GeV fixed, and considered two values 
of tan/? (0.5 and 1.0), and a range of values of M%. The resulting angles a and a& are 
rather independent of Mi as well as the choice of tan (3, whereas a c has some dependence 
on tan f3, as shown in the right panel of Fig. 7. 

For a given value of M 1; the resulting maximum is close to that found in [5], maximizing 
only with respect to the Hi contribution. We note that, considered as a function of M 1? 
there is a peak associated with the ti threshold. This is due to the contribution of the ti 
triangle diagram [5,7]. 

As discussed in [5], the heavier Higgs states have a tendency to reduce the CP- violating 
effect of the lightest one, unless they are sufficiently heavy to decouple. Thus, for a fixed 
value of the lightest Higgs mass, Mi, the over-all CP-nonconservation should increase as 
the second Higgs boson becomes heavier. This effect is illustrated in Fig. 8 for the case 
of Mi = 100 GeV and two values of tan/3 (0.5 and 1.0). Apart from some wiggles due to 
numerical noise, it is seen that there is a rather smooth increase of the sensitivity as the 
mass gap M 2 — Mi increases. 
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Figure 8: Maximum sensitivity S/N [see (5.5)] vs. M 2 , for fixed Mi 
parameters: M H ± = 700 GeV, £ per t = 5. 
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Let us now comment on the maximum CP nonconservation in the Yukawa sector, as 
given by the sensitivity in the quantity A±, compared with that of the gauge-Higgs sector, 
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£y. We already stated that these concepts are different. This statement can be made 
quantitative by considering the value of £y that corresponds to the rotation angles a, ab 
and a c for which the sensitivity in A\ is maximal. We find that £y ~ 0.6 and 0.3, for 
tan/5 = 0.5 and 1.0, respectively. 

6 Concluding remarks 

The concept of maximal CP nonconservation has been extended from the gauge-Higgs 
sector to the Yukawa sector, where various measures for CP nonconservation have been 
introduced and investigated. Large values of a and \a\,\ ~ 7r/4 favour large CP nonconser- 
vation in the Yukawa sector. But, in general, the maxima of CP nonconservation will in 
these two sectors not coincide. There could even be maximal CP nonconservation in one 
sector, and little or none in the other. 

We have here studied the simplest version of the 2HDM that allows for CP noncon- 
servation, where this CP nonconservation is given by one parameter, namely Im A5 in the 
potential (2.1). One could consider two more, independent parameters in the Higgs poten- 
tial that generate CP nonconservation, namely Im A 6 and ImA 7 (see, e.g., [6]). These terms 
in the potential are often considered less attractive, since they violate the Z 2 symmetry of 
the potential by terms which are quartic in the Higgs fields and thus make it more difficult 
to control flavour-changing neutral currents [13, 14]. 

However, if present, such terms would lead to a less constrained theory. While the 
Yukawa couplings (for Model II) are still given by the same elements of the rotation matrix 
R (and hence by the same expression in terms of tan (3 and the rotation angles a, ab and 
a c ), the masses M 2 and M 3 would be less constrained. By making these masses larger, the 
contribution of the lightest one, Hi, would be a better approximation to the over-all CP 
nonconservat ion . 

Acknowledgments: It is a pleasure to thank the organizers of the Epiphany 2003 Con- 
ference for creating a most stimulating atmosphere, and for excellent hospitality. 
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Appendix A. Maximizing 7^ 

This appendix deals with the maximization of 7y, Eq. (4.13). We shall first rewrite 7 & 7 t 
in terms of double angles. Let 

X = S b C c ~\- 5 C ) S b C c 5 C ) , 

V = (C 5 S b S c + S 5 C C )(S 5 S b S c - C & C C ), (A.l) 

then 

%% = z 2 (A.2) 

with 

z = 4c & s a s 2 b (c c s c ) 2 xy. (A.3) 

Maximizing amounts to maximizing the absolute value of z. 
We first note that 

x = |s 2 q[(1 + s 2 b )c 2c - c 2 b ] + \c 26l s b s 2c (A.4) 

where c 2& = cos(2d), c 2c = cos(2o; c ), etc. Furthermore, y can be obtained from x by the 
substitutions c & <-> Sq, and c c <-> s c , implying c 2(5 , <-> — c 2(5 ,, c 2c <-> — c 2c , with s 2(5 , and s 2c 
unchanged. Thus, 

Xy = {-\S2& C 2 b + [\c 2 a S b S 2c + jS 2& (l + Sfe) C 2c]} 
X { — \ s 2a C 2 — [\c 2 a s b S 2c + |s 2 q(1 + Sft)c2 c ]} 

= ]^[4 C b - 4 4 S b 4 - S L(! + S fc) 2c 2c - 4c 2a S 2 fi(l + sg) S b C 2c S 2c ] (A.5) 

The maximum is given by the three conditions: 

£ = 0, £- = 0, |1 = 0, (A.6) 
oa oath oa r 



or equivalently: 



3c 2fi (1 - 4)(1 - 4)(1 + s£) +4s 25 c 2c s 2c s b (l - 34) (1 + sg) 
+ 2c 2fi s 2 b [l - 74 - 34(1 - 34)] = 0, (A.7) 
(1 - 4)(1 - 4)(1 - 6s 6 b ) + c 2& s 2& c 2c s 2c s b (22st + 8s 2 b - 6) - 8^(1 - 44) 
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+ 413 - 27(44 + 7(4, + cl)} = 0, (A.8) 

C 2c (1 - <4)(1 - c4)(l + Sj) + C 2 a S 2 a «2c S 6 (l - 4<4)(1 + S*) 

-2c 2c ^[2c^ + 4(l-3c^)]=0 (A.9) 

While these three equations are highly non-linear, the solution of interest is actually 
obtained quite simply by setting 

C2a = 0, c 2c = 0, (A.10) 

whereby Eqs. (A. 7) and (A.9) become trivially satisfied, and Eq. (A.8) takes the simple 
form 

Qsl - 134 + - 1 = 0, (A.ll) 

the interesting solution of which is sf = 1/6. 
Summarizing, the maxima are obtained for 

a = ±±vr, a b = ±arcsiny/| = ±0.13386tt (24.1°), a c = ±\ir, (A.12) 

at which point 

7. = -V 

"8 x 1024 x 27 2 

determines the 70 of (4.17). 
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